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PREFACE 


As  part  of  The  Rand  Corporation’s  continuing  efforts  for  the  Ad¬ 
vanced  Research  Projects  Agency  In  the  study  of  laser  propagation,  a 
handbook  Is  being  prepared  that  presents  data  and  methodology  required 
for  calculating  the  atmospheric  degradation  of  laser  systems  in  tac¬ 
tical  missions.  The  theory  developed  in  this  report  will  be  an  essen¬ 
tial  Input  to  the  preparation  of  this  laser  propagation  handbook. 

The  outputs  of  many  tactical  lasers  are  best  approximated  by  a 
focused,  truncated-gausslan  disturbance  across  a  circular  aperture. 

A  quantitative  understanding  of  the  manner  in  which  an  initially  co¬ 
herent  beam  of  finite  cross  section  propagates  through  a  turbulent  at¬ 
mosphere  is  required  for  prediction  of  the  performance  of  various  de¬ 
vices  employing  lasers.  This  report  calculates  the  mean  intensity  dis¬ 
tribution  for  a  focused  gausslan  disturbance  in  a  finite  aperture  in 
both  the  near  and  far  field  and  examines  in  detail  the  effects  of  trun¬ 
cation,  wavelength,  and  aperture  size. 

These  results  should  be  of  use  to  those  Interested  in  propagation 
theory  and  its  applications  to  laser  range  finders,  laser  line  scan¬ 
ners,  communication  systems,  and  various  guidance  and  other  systems 
employing  an  illuminating  beam. 


SUMMARY 


The  first  section  of  this  report  briefly  outlines  the  development 
of  RM-6055-ARPA,  Propagation  of  a  Finite  Optical  Beam  in  an  Inhomoge¬ 
neous  Medium.  A  formula  is  derived  for  the  mean  intensity  distribu¬ 
tion  from  a  finite  beam  in  terms  of  complex  disturbance  in  the 
aperture  and  the  mutual  coherence  function  (MCFj  j.^.  w'"Heal  wave 
in  the  medium.  The  formula  is  used  to  examine  the  effects  of  turbu¬ 
lence  on  the  long-term  average  intensity  produced  by  a  focused,  trun¬ 
cated  gaussian  aperture  distribution.  In  particular,  it  is  shown  that 
(a)  while  the  vacuum  focal  point  intensity  will  increase  as  the  degree 
of  truncation  decreases  for  a  given  laser  output  power,  the  effect  of 
turbulence  limits  this  increase,  and  (b)  the  turbulence  can  virtually 
eliminate  the  vacuum  advantage  of  visible  over  infrared  wavelengths 
in  focusing  the  beam  at  practical  ranges.  Transverse  beam  patterns 
and  the  on-axis  intensity  are  shown  .for  CC^  wavelength,  and  a  crite¬ 
rion  is  established  for  the  condition  under  which  the  turbulence  pre¬ 
vents  effective  focusing. 
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I.  INTRODUCTION 


The  output  of  a  laser,  with  a  mirror  system  designed  to  keep  dif¬ 
fraction  losses  small,  can  usually  be  approximated  by  a  gaussian  ampli¬ 
tude  and  uniform  phase  distribution^  across  a  circular  transmitting 
aperture. *  Generally,  in  order  to  deliver  maximum  intensity  to  a  spec¬ 
ified  receiver,  the  output  beam  will  be  focused  onto  the  target,  with 
the  focusing  optics  truncating  the  gaussian.  The  complex  field  at  the 
exit  pupil  of  the  transmitter  for  this  configuration  may  be  approxi¬ 
mated  by  a  truncated-gaussian  amplitude  distribution  and  quadratic 
phase  distribution.  In  this  report,  we  will  calculate  the  mean  inten¬ 
sity  distribution  in  space,  resulting  from  the  above  complex  field  in 
a  transmitting  aperture  in  the  presence  of  a  turbulent  atmosphere. 

Previous  analyses  for  truncated  distributions  have  been  limited 
(2  3) 

to  vacuum  calculations.  *  The  atmospheric  calculations  have  been 
treated^'^’*^  by  introducing  the  fluctuation  in  refractive  index  as 
a  small  parameter  in  the  wave  equation,  and  expanding  the  fields  in 
powers  of  this  same  parameter.  This  procedure  results  in  the  field  be¬ 
ing  expressed  as  a  hierarchy  of  Born-type  integrals  over  the  vacuum 
fields  existing  at  each  point  in  space.  Because  of  the  extreme  com¬ 
plexity  of  these  integrals,  the  atmospheric  calculations  have  been 
limited  to  non-truncated  (possibly  focused)  gaussians,  for  which  a 
closed-form  solution  exists  for  the  vacuum  fields.  Even  in  this  spe¬ 
cial  case,  because  of  the  approximations  necessary  to  manipulate  the 
expressions,  it  is  not  difficult  to  show  that  none  of  the  expressions 
contained  in  Refs.  4,  5,  and  6  conserve  the  long-term  average  irradi- 
ance  over  a  plane  normal  to  the  direction  of  propagation,  implying 
that  even  those  solutions  must  be  incorrect. 

The  present  analysis  is  based  on  a  proof  developed  in  an  earlier 

(7)  ** 

study,  '  which  states  that  the  Huygens -Fresnel  principle  can  be  ex¬ 
tended  to  a  medium  exhibiting  a  spatial  variation  in  refractive  index. 
From  this  principle,  the  field  due  to  a  disturbance  specified  over  an 


With  confocal  mirrors,  the  lowest-order  mode  is  exactly  gaussian. 
See  p.  370  of  Ref.  8. 
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aperture  can  be  computed  by  superimposing  spherical  wavelets  radiating 
from  all  elements  of  the  aperture.  The  mean  intensity  distribution  is 
then  found  by  first  computing  the  intensity  at  a  point  from  an  arbi¬ 
trary  pair  of  elements  in  the  aperture.  Applying  a  reciprocity  theorem 
(also  proven  in  Ref.  7),  and  averaging,  reveals  that  the  above  quantity 
is  essentially  the  mutual  coherence  function  (or  MCF)  for  a  spherical 
wave  in  the  medium.  The  integration  over  the  aperture  is  performed  as 
a  final  step,  resulting  in  a  simple  formula  for  the  mean  intensity  pat¬ 
tern  valid  in  both  the  Fresnel  and  Fraunhofer  regions  of  the  aperture. 
The  development  is  outlined  in  Section  II. 

In  contrast  with  previous  work,  the  properties  of  the  medium  (e.g., 
the  turbulence  parameters)  appear  only  in  the  MCF  of  a  spherical  wave. 
The  present  calculation  separates  the  geometry  of  the  problem  (l.e., 
the  complex  aperture  distribution)  from  the  propagation  problem,  which 
is  determined  by  the  manner  in  which  a  spherical  wave  propagates  in  the 
medium.  The  MCF  describes  the  reduction  in  lateral  coherence  between 
different  elements  of  the  aperture,  effectively  transforming  it  into  a 
partially  coherent  radiator,  with  the  degree  of  coherence  decreasing 
as  the  distance  from  the  aperture  increases.  In  Section  III,  the  gen¬ 
eral  formula  relating  the  MCF  to  the  spectrum  of  index  of  refraction 
fluctuations  for  homogeneous  isotropic  turbulence  is  presented,  and  ap¬ 
proximate  formulas  are  given  for  the  MCF  at  various  propagation  ranges 
for  a  modified  Kolmogorov  spectrum.  A  curve  of  normalized  coherence 

length,  p  ,  defined  as  the  transverse  separation  at  which  the  coherence 
°  -1 

is  reduced  to  e  ,  versus  normalized  range,  is  also  presented. 

In  Section  IV,  the  case  of  a  focused,  truncated-gausslan  field 
distribution  is  examined  in  detail.  The  general  features  of  the  vacuum 
pattern,  with  particular  emphasis  on  the  effect  of  truncation,  are  dis¬ 
cussed  first.  For  a  fixed  laser  output  power,  a  curve  is  presented 
showing  the  vacuum  variation  in  focal-point  intensity  as  a  function  of 
the  degree  of  truncation.  The  curve  shows,  for  example,  that  when  the 

intensity  at  the  circumference  of  the  transmitting  aperture  is  reduced 
-3 

to  e  from  the  intensity  at  the  center,  the  intensity  at  the  focus  is 
reduced  to  e*85  percent  of  that  corresponding  to  a  uniform  distribution 
(this  latter  distribution  maximizes  the  vacuum  gain) .  The  effect  of 


turbulence  Is  shown  to  reduce  this  difference  in  all  cases  and,  in 
fact,  the  distribution  that  maximizes  the  gain  is  no  longer  a  uniform 
one.  A  criterion  is  also  given  for  which  no  further  increase  in  the 
on-axis  intensity  at  a  given  range  can  be  realized  through  focusing 
because  of  the  turbulence.  In  a  further  example  it  is  shown  that  for 
a  30-cm  beam  focused  at  0.5  km,  the  effect  of  the  turbulence  is  vir¬ 
tually  to  eliminate  the  vacuum  advantage  of  visible  over  infrared 
wavelengths  in  focusing  the  beam. 

In  addition,  beam  patterns  and  the  on-axis  intensity  are  shown 
for  CO2  wavelength,  for  both  a  30-cm  and  1-m  aperture,  using  turbu¬ 
lence  parameters  that  are  characteristic  of  paths  of  the  order  of  a 
few  meters  above  the  ground. 

The  results  are  briefly  summarized  in  Section  V,  where  it  is  ar¬ 
gued  that  the  spherical-wave  MCF  is  the  basic  quantity  to  be  measured 
for  computing  the  degradation  of  an  Intensity  pattern. 
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II.  THE  HUYGENS  FORMULATION 


Reference  7  shows  that  for  optical  propagation  In  a  weakly  In¬ 
homogeneous,  non-absorbing  medium,  the  field  at  the  observation  point 
P  from  a  complex  aperture  disturbance  U^(r^)  can  be  written  as 

w  ■  ir  J  G^*Ei)uA(£i)  d2si  (1) 

where  k  ■  2ir /X  Is  the  wave  number,  r^  Is  a  coordinate  In  the  plane  of 
the  aperture,  and  the  Integration  Is  carried  out  over  the  aperture. 

The  term  G(P,r^)  Is  the  field  at  point  P  due  to  the  spherical  wave 
propagating  In  the  medium  from  a  delta  function  source  at  r^.  Equa¬ 
tion  (1)  Is  valid  for  atmospheric  optical  propagation  where  the  scat¬ 
tered  field  can  be  shown  to  vary  slowly  over  a  wavelength  for  all 
propagation  distances  of  Interest  (e.g.,  wlO®  km  over  horizontal 
paths), ^  and  for  sufficiently  small  scattering  angles  where  the 
obliquity  factor'  '  can  be  taken  as  -1/X. 

The  Intensity  at  the  point  P  Is  then  given  by 

I(P)  -  U(P)U*(P)  -(^:)\(G(P,r1)G*(P,r2)UA(r1)U*(r2)  d2^  d2r2 

(2) 

It  Is  shown  In  Ref.  7  that  In  a  refractive  medium,*  the  complex 
field  at  P  due  to  a  point  source  at  r^  Is  Identically  equal  to  the 
field  at  r.  due  to  a  point  source  at  P— l.e.,  reciprocity  exists  In 
the  form  G(P,r^)  -  GCr^P).  Hence,  G(P,r^)G  (P,r2)  -  G(r^,P)G  (r2,P) 
Is  the  field  at  r^,  multiplied  by  the  complex  conjugates  of  the  field 
at  r2,  due  to  a  spherical-wave  source  at  P.  The  ensemble  average  (de¬ 
noted  by  angular  brackets) 


*[Iei-*I -I  £2-?l] 

<0(t.,?)C  Cr,,P)>  5  - - — - -  M(r,r,,P)  (3) 

I*!  -  El  ll2  -  El 


A  more  general  theorem  Is  proven  by  Rayleigh  (see  p.  380  of  Ref.  9). 


is  the  cross  correlation  of  the  complex  fields  at  the  points  r^,!^  ^ue 
to  a  unit  point  source  at  P.  When  the  scattering  angles  are  suffi¬ 
ciently  small  so  that  the  aperture  points  £^*£2  can  be  considered  as 
lying  on  the  surface  on  a  sphere  centered  at  P,  then  the  function 
MgCr^i^'P)*  defined  by  factoring  out  the  vacuum  fields  in  Eq.  (3),  is 
the  mutual  coherence  function  for  a  spherical  wave. 

If  the  vector  g  is  defined  as  the  normal  from  the  z-axis  of  sym¬ 
metry  to  the  observation  point  P,  then  in  the  small-angle  approximation 


r  (*i  -  2) 2 

e*l£i  -  elk[Z  ^  — 5r- 

Ie±  -  p!  “ 


(4) 


Hence,  we  obtain  from  Eq.  (2) 


"  (zT?)  j-  •  (El  *  E2>  +  rl  -  r2]j 


«  VEi-E2.J>»A<Ei>uI<E2>  ^Ei  «2E2 


(5) 


Changing  variables  in  Eq.  (5)  to  £  ■  (r^  -  rj) ,  r  *  £(r^  +  £2) » 
and  restricting  the  analysis  to  the  case  of  homogeneous  isotropic  tur¬ 
bulence  where  Mg^.rg.z)  ■  Mg(p,z)  (with  p  -  |g|),  yield 


(I>(g,z)  -  (£)2J  d2D  Ms(o,«)e"(lk/z)E’E  J«A(r  +  h_) 


.  U*(t  -  *cJe<lk/z>»’E  d2r  (6) 

Equation  (6)  gives  the  mean  spatial  intensity  distribution  in  both 
the  Fresnel  and  Fraunhofer  regions  from  an  arbitrary  complex  disturbance 
in  a  finite  aperture  in  the  presence  of  a  homogeneous,  isotropic,  tur¬ 
bulent  medium.  It  should  be  noted  that  integration  of  Eq.  (6)  over  the 

2  2  2 

g  plane  yields  /  <I)  d  g  •  f  |UA(r)|z  dzr.  Hence,  this  development  im¬ 
plies  energy  conservation  for  arbitrary  aperture  functions. 


/ 
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III.  THE  SPHERICAL-WAVE  MCF 


The  long-term  effect  of  the  turbulence  on  the  beam  pattern  enters 
through  the  spherical-wave  MCF.  Over  a  uniformly  turbulent  path.  Mg 
is  related  to  the  spectrum  of  index  of  refraction  fluctuations  by  the 
formula (10) 


Mg(p,z) 


dKK*n(K)  J*j  Jq(Kpu)  du" 


*n(K)K  dK 


(7) 


is  the  propagation  distance  at  which  the  mean  spherical-wave  field  is 
reduced  to  e  *  of  the  vacuum  field. 

Lutomirskl  and  Yura^1^  have  used  a  modified  Kolmogorov  spectrum 
for  the  turbulence 


yio 


2  -CKO 

0.033CS 
(K2  +  I'2) 11/6 


(9) 


where  C  is  the  index  structure  constant,  and  l  ■  2ir t  and  L  «  2irl 
n  o  o  o  o 

are  the  inner  and  outer  scales  of  turbulence,  respectively. 

In  Fig.  1,  a  plot  of  z  versus  X  is  shown  for  three  values  of 

2  c  .u  -2/3 

CR  roughly  corresponding  to  weak  (3  x  10  cm  '  ),  medium  (3  x 

10"15  cm"2^),  and  strong  (3  *  10”**  cm”2^)  turbulence.  For  hori¬ 
zontal  propagation  near  the  ground,  we  have  used  the  nominal  values 

t  ■  0.1  cm  and  l  ■  100  cm.  In  the  usual  atmospheric  case,  l  « 
o  o  o 

Lq,  and  the  integral  in  Eq.  (8)  can  be  approximated  to  yield 
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z  -  (0.39k2cV/3)  1  (10) 

c  no 

The  mean  field-decay  length,  z£,  and  hence  the  MCF,  will  depend  strongly 
on  the  outer  scale  of  turbulence. 

Reference  10  shows  that  there  are  three  distinct  propagation  ranges 

for  which  approximate  expressions  exist  for  the  MCF. 

The  first  is  found  by  observing  that  for  any  z,  M.  is  a  monotonl- 

*  -2z/z 

cally  decreasing  function  of  p  with  Mg(0,z)  -  1  and  Mg(«*,z)  -  e  c. 
Substituting  Eq.  (9)  into  Eq.  (7),  it  can  be  seen  that  this  asymptotic 
value  is  reached  when  p  »  Lq.  Physically,  when  the  separation  is  large 
compared  with  the  distance  over  which  the  refractive  index  fluctuations 
(or  temperature  fluctuations)  are  correlated,  ~Lq ,  then 

-2z/z 

Ms(p,z)  -  <U(r1,*)U(r1  +  p,z)>  -►  (U^,*))^^  +  £,z)>  -  e  c 

p  »  L 
o 


because  the  light  reaching  the  points  (r^.z),  (r^  +  g,z)  has  propagated 
through  essentially  statistically  Independent  media.  Hence  we  define 
Mg1*  ■  Mg(j>,  *  «  *c)  ■  1. 

To  describe  the  remaining  two  regions,  it  is  convenient  to  first 

characterize  the  MCF  at  a  given  range  by  the  separation  pq,  for  which 

Mg(pQ,z)  ■  e  1.  Then  the  second  region  is  found  by  observing  that  for 

l  «  p  «  L  ,  the  Integrand  in  Eq.  (7)  can  be  expanded  in  powers  of 
o  o  o 

(p/iQ)  to  yield'*' 

1 

■ |-1-‘(rXt)5/3[l  -  °-71(|-)1/3]|  <ll) 

(  co  L  0  J 

The  expansion  proves  valid  only  when  Mg(lQ,z)  «  1  and  >^(Lo,z)  «  1, 
which  implies  the  range  limitation 


(p,z)  -  exp  <-0.55k  Cnzp‘ 


i  -  0.7l(|.) 
L  o 
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z  «  z  «  z 
c 


where 


z 


i 


(0.39k2cV/3)_1 

n  o 


I 


For  ranges  greater  than  z^,  Mg(lQ,z)  <<  1,  and  all  of  the  p's  of  in¬ 
terest  are  small  compared  with  the  inner  scale.  Then,  for  pq  «  lQ, 
the  Bessel  function  in  Eq.  (7)  can  be  expanded  in  powers  of  (p/*0) 
to  yield 


M^(p,z)  -  exp  £-k2qzp2] 

[-•"  w’] 


■  exp 


z  »  zt 


(12) 


where 


q  -  <w2/3)  \  ♦  (K)K3  dK  w  0.56C2i-1/3 
n  no 


Hence,  for  propagation  paths  that  are  short  compared  with  z^,  the 
MCF  does  not  depend  on  the  inner  scale,  and  can  be  written  as 


o  c 


1  5 

For  the  range  0  <  z  «  z.  (for  L  fl  ■  10  ,  z.  -  10  z  ),  p  /fc  versus 

X  OO  1  coo 

z/z  ,  obtained  bv  inverting  the  equation  F(p  /£  ,  z/z  )  -  e-*,  is 
c  o  o  c 

plotted  in  Fig.  2.  For  ranges  z  <  £z  ,  M_(»,z)  >  e_1,  and  the  coher- 

C  j 

ence  length  as  defined  here  is  infinite. 


'I 


It  is  apparent  that  for  pQ  much  greater  than  the  aperture  diameter 
D,  Mg  in  Eq.  (6)  can  be  replaced  by  unity,  and  the  mean  intensity  is 
given  by  its  vacuum  value.  When  pQ  <  D,  the  turbulence  reduces  the 
average  irradiance  on  a  target. 


ii  nl  i 


i  i  i  m  i  hi  i  i  i  i  1 1  nl  i  i  i  i 1 1 1 


Fig.  2— -The  normalized  coherence  length,  p  fit  >  a*  a  function 

O  o 

of  z/z  for  a  spherical  wave 
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IV.  THE  FOCUSED.  TRUNCATED-GAUSS IAN  DISTURBANCE 

THE  MEAN  INTENSITY 

The  complex  field  in  the  aperture  corresponding  to  a  gaussian 
amplitude  distribution  with  standard  deviation  a,  focused  at  range  f, 
and  normalized  to  unit  amplitude  at  the  center  of  the  aperture  is 


"a  <hil> 


(a  ^+ikf  ^ ) 
e 


\r±\  €  D/2 


-  0 


Ie±I  >  d/2 


(14) 


With  given  by  Eq.  (14),  the  inner  integral  in  Eq.  (6)  is  the 
integration  of  the  function 


exp  [-a  2(r2  +  £p2)  +  iKz-^  -  f  1)p  *  r] 


over  the  area  of  overlap  of  two  circles,  each  of  diameter  D,  with  cen¬ 
ters  relatively  displaced  at  opposite  ends  of  the  vector  p.  The  inte- 

2 

gration  is  straightforward  and  yields  for  the  inner  integral  D  T.  „(x), 

0  ,  p 

where 


■  * 


-26  x 


2„2  />cos_1(x)  fl 

)o  *U_ 

COSCp 


2  2 

-6  (p  -2yx  cos<p) 


x  cos  [2Bx(u  cos  (p  —  x) ] u  du  x  <  1 


-  0 


X  >  1 


(15) 


2  —1  —1 

where  x  ■  p/D,  6  ■  D/2a,  and  (5  ■  ^kD  (z  -  f  ) .  Then,  using  polar 
coordinates  for  the  p  integration,  performing  the  angular  integral, 
and  changing  variables  from  p  to  x  ■  p/D,  yield 
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(l)(a,z)  -  f  (^")  ^  5cJo(2ax)Ms(Dx,z)r6tB(x)  dx 


where  JQ  Is  the  zero-order  Bessel  function,  and  a  ■  kDp/2z  ■  fkD  tan  8, 
where  6  Is  the  angle  that  the  direction  to  P  makes  with  the  central  di¬ 
rection. 

THE  VACUUM  PATTERN  AND  EFFECT  OF  TRUNCATION 

In  the  limit  of  zero  turbulence,  Mg  ■  1,  and  Eq.  (16)  reproduces 
the  vacuum  Intensity  pattern.  In  this  limit,  it  follows  from  the  def¬ 
inition  of  0  that  the  angular  pattern  In  the  focal  plane  (z  ■  f)  is 
identical  with  the  Fraunhofer  pattern  for  the  unfocused  case  (f  -  “, 
z  ■  ») .  Further,  because  r.  a(x)  Is  an  even  function  of  3,  the  angular 
vacuum  patterns  are  Identical  In  planes  for  which  |z  1  -  f  x|  -  constant. 

For  the  limit  D  «  a  (6  -*  0) ,  the  problem  reduces  to  this  focusing 
of  a  plane  wave  diffracted  by  a  circular  aperture.  No  simple  expres¬ 
sion  exists  for  the  transverse  vacuum  Intensity  except  In  the  focal 
plane,  where  the  distribution  is  given  by  the  Airy  pattern 


,»•■>•  t(#)M 


The  on-axls  Intensity  is  given  by 


„  .  1  /kD2\2/Sta  4^ 

,<°.o  •  7  \xr)  ^~]f~ j 


which  has  an  absolute  maximum  at  the  focal  point.  The  on-axls  nulls 
are  located  at 


16nirf 


n  ■  1,2, . . . 
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2 

When  16irf/kD  >  1,  Che  only  on-axis  nulls  lie  beCween  the  transmitter 

2 

and  the  focal  plane;  when  16irf/kD  <  1,  there  are  on-axis  nulls  beyond 

2 

the  focal  plane.  In  the  limit  16irf/kD  «  1,  the  nulls— indeed  the 
whole  pattern— are  symmetric  about  the  focal  plane  for  ranges 
|(z  -  f)/f|  «  1. 

In  the  limit  D  »  a  (or  {-►<»),  when  the  aperture  does  not  effec¬ 
tively  truncate  the  gaussian(  the  vacuum  pattern  propagates  as  a  gaus- 
sian  in  all  transverse  planes  according  to  the  equation 

i  2.2  2 

I  (p,z)  -  ^  e‘P  /a  W  (19) 

8  w 

where 


The  maximum  intensity  in  this  case  does  not  occur  at  the  focal  point, 

2 

but  rather  on  the  axis  at  the  shorter  range  z  -  f/[l  +  (f/ka  )],  where 

2  2 

the  intensity  is  1  +  (ka  / f)  .  Generally,  in  order  to  focus  the  laser 

beam  effectively,  the  focal  length  should  be  much  less  than  the  smaller 
2  2 

of  ka  ,  kD  ,  in  which  case  the  intensity  at  the  focal  point  is  approx¬ 
imately  the  smaller  of  (ka2/f)2,  ^(kD2/f)2. 

For  the  case  of  arbitrary  truncation,  the  on-axis  intensity  can 
be  shown  to  be 


with  no  simple  formula  for  the  transverse  pattern,  even  in  the  focal 
plane.  The  average  aperture  Intensity,  defined  as  the  total  laser  out¬ 
put  power,  W,  divided  by  the  area  of  the  aperture,  can  be  computed  by 
Integrating  | U^J 2  over  the  aperture,  and  is 
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(21) 


The  effect  of  truncation  on  the  intensity  at  the  focal  point  for 
a  given  output  power  can  then  be  computed  from  the  gain 


It  has  been  shown^^^  that  of  all  possible  aperture  distributions, 
the  focal-point  intensity  for  a  given  W  is  maximized  by  a  uniform  dis¬ 
turbance,  corresponding  in  our  case  to  6  ■  0,  and  yielding 


Fig.  3 — The  normalized  gain 
as  a  function  of  82 

the  intensity  at  the  center 
vacuum  reduction  of  **15  percent 
shown  to  reduce  this  difference 


In  Fig.  3  the  normalized  gain 


2 

is  plotted  versus  5  ,  the  power 
of  e  by  which  the  intensity  at 
the  circumference  is  reduced  from 


.3 

For  example,  a  ratio  of  e  yields  a 
The  effect  of  turbulence  will  be 


THE  EFFECTS  OF  TURBULENCE 

The  long-term  effect  of  the  turbulence,  as  discussed  in  Section 
III,  is  to  reduce  the  lateral  coherence  between  the  fields  radiating 
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from  different  elements  of  the  aperture  for  separations  pq.  For 

pQ  <  a,D,  the  beam  will  be  confined  to  a  circularly  symmetric  area 

of  dimension  ~z/kp_  at  a  distance  z.  The  ratio  of  the  field  intensity 
0  2 

to  aperture  intensity  is  then  ~(kDpQ/z)  ,  and  the  on-axls  intensity 
can  be  Increased  by  focusing  the  beam  only  if  this  quantitly  is  con¬ 
siderably  greater  than  unity  at  the  range  of  interest.  In  general, 
this  partially  coherent  aperture  will  produce  an  intensity  pattern 
with  less  detailed  structure  than  the  vacuum  distribution,  and  with 

the  maxima  and  minima  tending  to  be  averaged  out.  Further,  because 
2 

z  %  «  X  ,  there  will  be  less  degradation  of  longer  wavelength  patterns. 

In  Figs.  4a-4c,  the  focal-plane  atmospheric  patterns,  normalized  to 
the  vacuum  focal-point  intensity  (1/4) (kD  /4f)  ,  are  shown  for  a  uni¬ 
form  aperture  disturbance  with  D/*>0  ■  0.01,  0.1,  and  1,  respectively. 
With  Z£  determined  from  Eq.  (10)  (or  Fig.  1),  these  curves  can  be  used 
for  arbitrary  wavelengths.  The  dashed  curve  in  each  figure  is  the 
vacuum  Airy  pattern,  Eq.  (17),  similarly  normalized.  For  a  given  ra¬ 
tio  of  z/z  ,  the  greater  the  diameter  (D /L  ),  the  less  the  effective 
c  o 

coherence  over  the  aperture,  and  hence,  the  greater  the  degradation 
from  the  vacuum  pattern. 

In  the  remaining  examples,  the  normalized  intensities  (I)N»  de¬ 
fined  as  the  ratio  of  the  mean  intensities  (I)  to  the  mean  aperture  in- 
_  2 

tensity  I  ■  4W/irD  ,  are  plotted.  Turbulence  parameters  have  been 

chosen  that  are  characteristic  of  the  strengths  and  scales  found  from 

2  -5  -2/3 

zero  to  several  meters  above  the  ground:  C  ■  3  x  10  cm  ,  !  ■ 

n  o 

0.1  cm,  and  ■  100  cm. 

Figures  5a  and  5b  display  the  effects  of  truncation  on  the  focal- 

plane  intensity  of  a  10.6-p  beam  from  a  1-m  aperture  focused  at  0.5  km. 

The  vacuum  and  atmospheric  patterns  are  shown  in  Fig.  5a  for  a  uniform 

distribution,  whereas  in  Fig.  5b  the  gausslan  Intensity  distribution 
-4 

is  reduced  to  e  at  the  circumference  from  its  value  at  the  center. 

In  vacuum  the  truncated  distribution  has  a  focal-point  Intensity  of 
«25  percent  less  than  the  uniform  case;  however,  the  effect  of  turbu¬ 
lence  is  to  reduce  both  intensities  such  that  the  difference  is  only 
«wl0  percent.  Because  the  MCF  reduces  the  effective  average  cohererce 
between  different  radiating  elements  of  the  aperture  in  a  manner  that 
depends  only  on  the  distance  between  them,  elements  near  the  edge  of  the 


Comparison  of  focal  plane  patterns  at  X  =  10.6  /x  for  D  =  1  m  focused  at  0.5  km 
for  (a)  uniform  and  (b)  truncated  gaussian  aperture  distribution 
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aperture  weigh  less  than  those  near  the  center  in  contributing  to  the 
average  Intensity.  Hence,  the  relative  intensity  gain  of  a  uniform 
distribution  over  a  truncated-gaussian  distribution  will  always  be 
less  than  in  the  vacuum  case.  In  fact,  it  can  be  shown  from  Eqs.  (6) 
and  (21)  that  in  the  presence  of  turbulence,  the  aperture  distribution 
which  maximizes  the  gain,  G  -  (l)(0,f)/I,  is,  in  general,  not  uniform, 
but  somewhat  peaked  at  the  center  of  the  aperture. 

In  Figs.  6a-6c,  the  focal-plane  vacuum  and  atmospheric  patterns 
are  compared  for  D  ■  30  cm,  6-0,  and  z  -  f  -  0.5  km,  using  the  wave¬ 
lengths  0.6328  y  (He-Ne) ,  1.06  y  (Nd) ,  and  10.6  y  (CO^) .  Although  the 

vacuum  intensities  at  the  focal  point  for  the  shorter  wavelengths  are 

_2 

greater  (the  "antenna  gain"  is  proportional  to  X  ) ,  the  effect  of  the 
atmosphere  for  the  chosen  geometry  essentially  eliminates  the  shorter 
wavelength  advantage  in  focusing  the  beam.  For  larger  apertures,  the 
advantage  of  CO2  in  producing  an  intense  focal  intensity  will  be  greater. 

In  Figs.  7a-7c,  for  X  -  10.6  y,  the  vacuum  and  atmospheric  focal- 
plane  patterns  are  shown  for  D  -  1  m  and  6-0  (uniform  distribution) 
at  focal  lengths  of  0.5  km,  5  km,  and  10  km.  The  long-term  effect  of 
the  turbulence  at  the  larger  ranges  is  to  reduce  the  on-axis  intensity, 
and  to  spread  out  the  beam.  Hence,  the  positioning  accuracy  required 
to  illuminate  a  target  can  be  considerably  less  than  a  vacuum  calcula¬ 
tion  would  indicate.  As  discussed  above,  the  solid  curves  (with  turbu¬ 
lence  present)  are  relatively  insensitive  to  the  degree  of  truncation. 

Finally,  in  Figs.  8a  and  8b,  the  vacuum  and  atmospheric  on-axis 
intensities  are  compared  for  a  10.6-y  beam  focused  at  0.5  km  with  6  ■ 

0,  and  D  -  30  cm  and  1  m,  respectively.  The  intensities  are  shown  for 
a  range  of  50  m  on  either  side  of  the  focal  plane.  The  on-axis  inten¬ 
sity  for  the  30-cm-aperture  case  remains  reasonably  constant  over  the 
100-m  range,  with  a  maximum  at  «30  m  on  the  aperture  side  of  the  focal 
plane.  The  axial  distribution  for  the  1-m  aperture  is  almost  symmetric 

with  respect  to  the  focal  plane  and  obviously  requires  greater  focusing 

4 

accuracy  to  produce  the  available  gain  of  wlO  at  a  specified  target. 

At  ranges  of  5  and  10  km,  the  on-axis  intensities  are  essentially  con¬ 
stant  within  50  m  of  the  focal  plane  and  are  equal  to  the  focal-point 
values  of  Figs.  7b  and  7c,  respectively. 
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V.  DISCUSSION  AND  CONCLUSIONS 


Based  on  the  results  of  Sections  III  and  IV,  the  following  pro¬ 
cedure  can  be  used  to  estimate  the  focal-point  Intensity. 

First,  from  a  knowledge  of  (1)  the  Index  structure  constant  Cr, 
(2)  the  (outer)  scale  over  which  the  temperature  fluctuations  are  rea¬ 
sonably  correlated,  Lq  and  (3)  the  wave-number  of  the  radiation,  k, 
the  mean-field  decay  length  can  be  constructed  according  to  Eq.  (10) 
as 


*  -  Uw'f 

c  \  no/ 

Then  the  equivalent  coherent  aperture  size,  p  ,  at  range  z,  can  be 

2 

computed  from  Fig.  2.  If  the  smaller  of  the  two  quantities  (kDp  /z)  , 

2  ° 
(kapQ/z)  is  considerably  greater  than  unity,  then  the  intensity  at 

the  range  z  can  be  increased  by  focusing  the  beam  at  that  range.  If 

this  Inequality  is  not  satisfied,  no  Intensity  amplification  can  be 

realized  by  focusing  the  beam. 

These  estimates  are  also  correct  for  the  case  of  Inhomogeneous 
turbulence  when  pQ,  corresponding  to  the  e  *  point  of  the  MCF  for  the 
Inhomogeneous  medium,  is  used.  The  formula  for  the  spherical-wave 
MCF  In  an  inhomogeneous  medium  is  given  in  Ref.  10.  In  either  case, 
if  a  more  precise  estimate  of  the  degradation  due  to  the  turbulence 
is  required,  the  Integral  of  Eq.  (16)  can  be  evaluated.  This  will 
usually  be  necessary  when  pQ  ~  a  or  D.  When  pq  «  a,D,  then  the  ap¬ 
proximate  formulas  of  Ref.  7  are  applicable. 

While  the  expressions  for  pQ  depend  upon  particular  models  for 
the  turbulence  spectrum  for  homogeneous  turbulence,  and  model  profiles 
for  the  Inhomogeneous  case,  the  intensity  distribution  can  always  be 
determined  from  Eq.  (6)  if  the  MCF  is  known.  Hence,  in  order  to  pre¬ 
dict  the  beam  pattern  from  an  arbitrary  wavefront  in  an  aperture,  it 
is  the  MCF  that  should  be  measured  rather  than  specific  beam  patterns. 
The  measurement  of  turbulence  parameters  that  rely  on  specific  models 
for  the  spectrum  are  also  of  limited  utility. 


-23- 


Thls  point  can  be  clarified  by  first  observing  that  Eq.  (16)  is 
the  Fourier-Bessel  transform  of  the  quantity 

(f)(!r)  Ms(Dx‘z)rs,6<x) 

which  can  be  inverted  to  yield 


d2i,«,$(d)  vp*z) " 2*  \0  j0(3sSE)<i><p  ,z)p  dn  (24) 

where  (I)(p,z)  is  the  mean  intensity  at  a  distance  p  from  the  axis  at 
range  z. 

For  an  arbitrary  disturbance  in  the  aperture,  the  function  ^(x) 
given  by  Eq.  (15)  would  be  replaced  by  the  appropriate  overlap  Integral 
of  Eq.  (6).  Equation  (24)  thus  provides  a  possible  method  for  deter¬ 
mining  part  of  the  MCF  from  measurements  of  the  beam  pattern. 

However,  because  g(p/D)  ■  0  for  p  >  D,  inverting  the  intensity 
distribution  can  give  no  Information  regarding  Mg(p,z)  for  spatial 
separations  larger  than  the  diameter  of  the  transmitting  aperture.  In 
particular,  in  order  to  determine  the  distance  z£  from  a  beam  pattern 
measurement,  it  would  be  necessary  to  have  an  aperture  diameter  greater 
than  the  largest  scale  of  turbulence,  Lq,  which  might  be  of  the  order 
of  meters.  Hence,  if  one  can  determine  the  spherical-wave  MCF  for  all 
spatial  separations  at  a  given  range  (e.g.«  by  using  an  Interferometer), 
one  can  infer  the  intensity  distribution  from  an  arbitrary  aperture 
distribution  at  that  range,  while  the  reverse  is  not  true  unless  aper¬ 
tures  greater  than  the  coherence  length  at  that  range  can  be  constructed. 
Even  if  the  beam  pattern  were  measured,  Mg(p,z)  would  first  have  to  be 
constructed  from  Eq.  (24)  from  the  given  measurement  to  determine  the 
general  response.  Therefore,  for  design  purposes,  direct  measurements 
of  the  MCF,  which  are  notably  lacking,  are  required. 
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